Three "quantum" models of competition and cooperation in interacting biological 

populations and social groups 
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In present paper we propose the consistent statistical approach which appropriate for a number of 
models describing both behavior of biological populations and various social groups interacting with 
each other. The approach proposed based on the ideas of quantum theory of open systems (QTOS) 
and allows one to account explicitly both discreteness of a system variables and their fluctuations 
near mean values. Therefore this approach can be applied also for the description of small populations 
where standard dynamical methods are failed. We study in detail three typical models of interaction 
between populations and groups: 1) antagonistic struggle between two populations 2) cooperation ( 
or, more precisely, obligatory mutualism) between two species 3) the formation of coalition between 
two feeble groups in their conflict with third one that is more powerful . The models considered in 
a sense are mutually complementary and include the most types of interaction between populations 
and groups. Besides this method can be generalized on the case of more complex models in statistical 
physics and also in ecology, sociology and other "soft' sciences. 

PACS numbers: 03.65.Ta 



I. INTRODUCTION 

Among various classical open systems of interest in 
physics and also in so-called "soft" sciences such as ecol- 
ogysociologyeconomics and so on, the important role is 
played by the systems whose states in accordance with 
the sense of the problem are specified by the set of inte- 
ger variables {rii} (where n, = 0, 1, 2, and i =1,2,. ...N 
, where N is a number of degrees of freedom). For example 
in statistical physics n, may represent occupation num- 
bers of various cells in phase space of the nonideal gas ,in 
ecology - numbers of individuals in distinct populations 
living in the area and interacting with each other, in eco- 
nomics the number of different companies operating on 
the market at the same time and so on. As far as states 
of such systems change with time they are described as a 
rule by a set of autonomous differential equations of the 
next form: 



drii 
~dt 



Fi({n a }), 



(1) 



where Fi({n a })— some nonlinear state depending func- 
tions determined by a concrete problem. 

Obviously, that notation of Eq. (Q]) assumes that all 
variables n, - are continuous quantities. In the case when 
all rii 1 such " smooth" approximation of discrete sys- 
tem is quite reasonable. But for small occupation numbers 
when m > 1, the dynamical approach of Eq. ([T]) becomes 
inapplicable and it is necessary to take into account both 
as discreteness of variables rii and their possible fluctua- 
tions near mean values. The statistical method proposed 
in the present paper has undoubted advantage compared 



with dynamic approach Eq. ([T]) since it completely free 
from mentioned restrictions. 

The method proposed is based on application of the 
quantum Lindblad master equation (LME) for density 
matrix (DM) evolution of a open quantum Markov sys- 
tem. The LME in general case has the next form( 



d P i \u 1 



■Y, [ftp, Rf]+h. 



(2) 
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where H is some hermitian operator, describing the inner 
dynamics of open system and {Ri} are a set of nonhermi- 
tian operators that simulate different types of interaction 
of open system in question with its environment). 

Although Eq. (|2J has obviously quantum origin, never- 
theless, under the certain conditions which we will discuss 
more detail later in this paper, it may be applied also to 
the study of many classical open systems ,for example for 
the systems with discrete variables. The main advantage 
of the LME compared with other similar master equa- 
tions lies in the fact that the number and form of opera- 
tors Ri entering in Eq. ([2]) for many concrete problems 
can be defined from simple heuristic reasons. As soon as 
such choice is made the consistent mathematical frame- 
work for the description of the problem is at our disposal. 
Since the method proposed is to some extent heuristic, we 
demonstrate its reasonableness and effectiveness on the 
examples of different typical models which equally can 
be applied to ecology for the description of behavior of 
interacting populations and in sociology for striding dy- 
namics of cooperating or competing groups. Although in 
the present paper we specially restricted our choice only 
the simplest examples, which admit the complete quali- 
tative analysis , it is clear that the method proposed can 
be generalized on more complex cases without any dif- 
ficulties of principle. The rest of the paper is organized 
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as follows. In Sect. 2 we present the main features of ap- 
proach proposed , and discuss all conditions and restric- 
tions for its application .In Sect. 3 we study the model of 
the struggle for existence between two antagonistic bio- 
logical populations in the case when total size of both 
populations remains invariable. In Sect. 4 we consider the 
model of cooperation (or more precisely the model of 
obligatory mutualism) between two distinct populations 
or social groups. In Sect. 5 we study somewhat more com- 
plicated model of three groups interaction in situation 
when two feeble groups join together to confront success- 
fully more powerful rival. In conclusion we are summing 
up all results obtained in the paper. 

Now let us turn to the concrete presentation of the 
paper. 



II. DESCRIPTION OF THE METHOD 

As we already note the mathematical framework of the 
method proposed is the Lindblad master equation ([2]) 
that is used in QTOS for the description of open quantum 
Markov systems. Although this equation has undoubtedly 
quantum origin under certain conditions it can be applied 
also for statistical description of some classical open sys- 
tems. 

This significant point should be explained more de- 
tail. First of all note that density matrix (DM) of quan- 
tum system apart from quantum correlations containes 
also exhaustive information about its classical correla- 
tions. We believe that information about classical corre- 
lations is recorded mainly in diagonal elements of DM. 
Let us assume now that the open system of interest is 
such that the relevant LME implies closed and complete 
set of equations connecting only diagonal elements of 
its DM . Clearly, in this case the LME could serve as 
appropriate framework for the statistical description of 
classical analog of corresponding quantum system. Al- 
though similar situation is realized for several classes 
of open classical systems, but in the present paper we 
are interested in only special such class, namely clas- 
sical systems, whose states can be specified by integer 
variables n*. Let us describe this class more explic- 
itly. Our first and main assumption is that " hamiltonian 
term" in the r.h.s. of the LME ([2]) is missing. Besides 
we demand that all operators Ri in Eq. @ have mono- 
mial form that is: Ri = d ]J ( a J) fe, '° (a/?)'* ' 3 ( m g en_ 

eral case Ci can be considered as functions of occupa- 
tion numbers , but in what follows they are assumed to 
be constants). In such situation one can directly verify 
that all equations for diagonal elements of /o(ni, ...njv), 
(where m — a^di and a,i,a~£ are usual Bose-operators, 
that is [aj,a^"] = Sik,) form complete set of differential- 
difference equations, which give self-consistent description 
of the system. It should be note also that with assump- 
tions made the LME ^ actually reduced to the specific 



form of Pauli master equation (PME) jH, namely: 

^ = Yl W »<*Pm ~ WnxnPn). (3) 

Tlx 

We could suggest also that for the" good enough" form 
of transition probabilities W nni entering in the Eq. Q 
one can find an appropriate set of operators Ri of mono- 
mial form such that the LME (J2J) (without "hamiltonian 
term") and the PME ^ actually will coincide. Now to 
complete the description of the method we need to spec- 
ify how one must choose a set of operators {Ri} for a 
concrete problem of interest. The quality of the method 
proposed in our opinion lies exactly in this point becouse 
this choice can be done according to simple heuristic con- 
siderations. Namely, we believe the number and form of a 
set {Ri} are entirely determined by the condition what 
types of transitions W nni one want to take into account 
for the concrete problem. The best way to illustrate all 
features of the method proposed is to apply it to study 
various concrete models .Let us now turn to this matter. 



III. ANTAGONISTIC STRUGGLE BETWEEN 
TWO COMPETING POPULATIONS 

As the first example that illustrates all features of the 
method proposed we consider the model of antagonis- 
tic opposition between two populations or social groups 
struggling with each other for certain resourses or some 
other preferences. In most naked form such struggle is 
realized when two tribes or kins of cannibals living side 
by side in the area in a literal sense eat each other ( but 
with unequal voracity). Assume for the simplicity that 
the total number of individuals in both populations re- 
mains constant in this struggle. As regards to sociology 
it could be for example the struggle between two polit- 
ical parties competing at the parlament elections when 
the total number of vacant seats is fixed. It is clear that 
in such situation the benefit for one group necessarily 
means the failure for the other and vice versa. Using the 
language of game theory, one can say that we consider 
statistical model of two person game with zero sum (but 
from nonstandard point of view). Let us turn now to 
the explicit mathematical formulation of the model. We 
assume that our model of antagonistic struggle can be 
properly described by the help of the LME with two op- 
erators: R\ = yj^a^a2 and R 2 = \J\ a 2 a i (where op- 
erator Ri corresponds to events when first population or 
group benefits and the second fails, and R 2 respectivly to 
the opposite one, coefficients a and b reflect correspond- 
ing competitiveless of both groups) . After these principal 
assumptions remainder of our analysis is completely rig- 
orous and can be represented as the sequence of three 
consecutive steps. 

Stepl. We write down the LME for the DM of the 
system , that characterize total correlations between two 
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groups of individuals with antagonistic interaction. The equation has the next form: 
I 



dp 
~db 



a (a^"a2/oa^ai — aia^ a2p) + 6 (a^~aipa^"a2 — a 2 pa\ ai) 

I 



(4) 



One can easily see that Eq.Q implies the closed set of only classical correlations of the model in question, 
equations for diagonal elements of DM that describes namely: 

I 



d P ni , 
dt 



n l n 2P ni n5 ~ n 2niP nin2 



n 2 mpj^ n2 - n in2 p nin2 



(5) 



(note that we use the notation n = n + 1, and n = n — 1 
to reduce the length of formula ([5]) 

Step2. It is convinient to represent the difference- 
differential system of equations (JSJ in the form of equiv- 
alent differential equation. Such representation can be 
obtained with the help of the generation function for 
the distribution p nin2 - According to definition, the gen- 
erating function (GF) G(u,v,t) = Pmn 2 univn2 ■ 

The normalization condition ^ P ni n 2 = hnplies that 



G(l,l,i) = 1. All moments of distribution P ni n 2 
can be found by differentiating G(u,v,t). For example: 

dG I 



ni 



= i,ni (ni - 1) 



d 2 G 



One can verify directly that Eq 
equivalent equation for GF G (u, v, t) 



| u= „=iand so on. 
(|3J) implies the next 



dG 
~dt 



a(u — v) 



d 2 (uG) 



b(v 



d 2 (vG) 



(6) 



dudv dudv 

Let us show that Eq. (J6)) admits an exhaustive qualitative 
analysis because its general solution can be represented 
in the form of the decomposition: 



G (u, v, t) — K N Gn (u, v, t) . 



N 



(7) 



where Gn (u>, v, t) = A^u + Aiu v + ...Ajyv is nor- 
malized homogeneous polynomial of degree N in vari- 
ables u and v; and {Kn} is a set of constants satisfying 
to the condition : = l.The capability of decom- 

JV 

position Eq. J7| is a consequence of two reasons: 1) the 
linearity of Eq. ([6]) and 2) the existence of conserva- 
tion law, namely, one can easily see that total number 
of individuals in both populations N = ni + ni- remains 
constant. 

Thus we conclude that for the study of general dy- 
namics of Eq. (jSJ) enough to consider it in each of the 
subspaces with fixed N, where this dynamics can be re- 
duced to the simple system of linear equations for the 



coefficients Ai of homogeneous polynomial G 



N ■ 



dAj 
dt 



LitAk- 



(8) 



In the Eq. © is some nonhermitian N + 1 x 
N + 1 matrix which can be easily calculated from Eq. 



For example in the sim- 
has the form: L 2 = 

Having in hands decompo- 



([6]) for any concrete N. 
plest case N = 2 matrix 
-26 2a 
2b - (2a + 26) 2a 
26 -2a J 

sition (|7|) we are able to answer on the main question: 
how the time evolution of arbitrary initial distribution 
Go (u, v) = Gn (u, v, t = 0) is happened? To answer this 
question let us begin with the end and find the stationary 
solutions of Eq. ([6]). Equating the r.h.s. of ((6]) to zero 
we obtain that: 



G st («,«) 



A (u) + B (v) 



au 



bv 



(9) 



To obtain required polynomial form of G (u, v, t) we must 
choose An (u) as Cn (au) +1 and B (v) as —Cn (bv) +1 
after that Eq. © leads to the next result: 



G# (u,v)=C N 



(au) 



N 



(au) N 1 (bv) 



..{buy 

(10) 

( where Cn = a N + a N 1 b + ...6 is normalization con- 
stant). With the help of Eq. flTDJl for the GF one can 
find all moments that is all statistical characteristics of 
the system of interest in its stationary states. In this con- 
nection we would like to mention one curious result. Let 
N is total ( and conserved) number of individual in both 
populations and assumed that a > 6. Then with the help 
of (fTT)|) one can show that when t tends to infinity, the 
fraction = tends to zero, that is less voracious popula- 
tion gradually disappears. 

This conclusion is completely coincides with result ob- 
tained from dynamical approach of Eq. (fTJ), since for 
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a 2 +ab+b 2 ' 

Going over to the 



N 1, the role of discreteness of variables n.j and 
their fluctuations becames negligible. On the other hand 
the expression (fTQ| is valid also in the case of small 
populations where dynamical approach is not applica- 
ble. For example when N — 2 GF is equal to: G2 = 

a2 ' U2 a2+ab+b* b2v2 EIld W(3 C£Ln filld th£Lt " 

^r— _ 2b 2 +ab ab 

" 2 — a 2 +ab+b 2 ' 111112 ~ a 2 +ab+b 2 - 

study of time evolution of arbitrary initial distribution 
Go {u, v) we propose the following precription: using Eq. 
one must expand Go in homogeneous polynomials 
Gjv, thus determining a set of Km , and after that we 
can immediately write down the definitive result, namely 
G(u,v, 00) = KnG n (u, v) .In what follows we will 

TV 

call states with GF Gn the pure states. 

Note also that using the Euler theorem about homoge- 
neous functions: u^f- + v ^jff~ = NGn , one can easily 
obtain for the pure states two simple relations: 



n\ + n\U2 = Nn\ and n 2 , + Ti\ni — Nn%. 



(11) 



If we define by standard way the variance of the num- 
ber of individuals in every population as Oi = nf — 
(ju) 2 (where i=l,2) and correlation coefficient k = 
Wnt-W ■ ^ then the Eq. (JTTJ implies that a x = <r 2 = 

— ri\ni and hence for any values a and b 
correlation coefficient k in any pure state is equal to 
— f .This result exactly justifies the name for considered 
model as antagonistic model of interaction between two 
populations or groups. 



tions either together get the benefit in the struggle for 
existence or under adverse circumstances are failing to- 
gether. In biology such kind of interaction is known as 
symbiosis ( or more precisely as obligatory mutualism). As 
the simplest example of such interaction we can specify 
the alliance between hermit- crab and actinium, that is 
sea coral. [4] Actinium attached to the shell in which crab 
lives , and then moves with it and eats the remaines of 
its food. On the other hand actinium protects the crab 
from its enemies with special cells located in its tenta- 
cles. Thus neither actinium neither hermit -crab cannot 
successfully exist without each other. As regards to soci- 
ology the number of distinct groups that share common 
goals and (or) values and cooperates with each other to 
achieve them are truly unlimited. To formulate the math- 
ematical model describing the case of two groups oblig- 
atory cooperation we start again from Lindblad equa- 
tion ([2]) , but now to accordance with the sence of the 
problem as a collection of operators Ri we choose two 

operators:i?i = and R2 = \J\a1a2 -The first 

operator describes the events in which, both groups get 
the benefit while the second one is responsible for the 
common failure. The required Lindblad equation for di- 
agonal elements of density matrix of the system p niTl2 in 
this case takes the form: 



IV. COOPERATION BETWEEN TWO 
POPULATIONS OR SOCIAL GROUPS 

In this part we consider the model of interaction be- 
tween two groups of individuals in which both popula- 



J 



dt 



nm 2 p ni na -ni-na- p nin2 



r 



■ b [m ■ 712 • {hfx ~ nW2P ni „ 2 ] 



(12) 



(In the Eq. (fT2j) as before we used the notation:?!! = 
m + 1 and rij = rij — 1). It is easy to verify that 
difference-differential equation (|12[) by introducing the 
generating function G (u, v, t) — p ni „,u Tll 7j™ 2 can be 

transformed to the equivalent form of single differential 
equation , namely: 



dG 
~dt 



(1 



d 2 



dudv 



[{b - auv) G] 



(13) 



Note that Eq. f|13|) as well as Eq. P^)l admits the exhaus- 
tive qualitative analysis, due to linearity and the presence 



of conservation law. We begin our analysis of Eq. (fT3")) 
with finding its stationary solutions. Equating the r.h.s. 
of Eq. (Ti~3"f to zero, we obtain that a set of stationary 
solutions of Eq. (|L3|) can be represented in the form: 



G st (u,v) 



A(u) + B(v) 



(14) 



where functions A (u) and B (v) yet to be determined. 
As the second step let us write down equations of motion 
for the occupation numbers mean values of both popu- 
lations. By differentiaiting of Eq. (|13p with respect to u 
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and v respectively and then putting u 



1 we obtain: 



dni 

at 



dn 2 

a (W[ + + 1) + (a - b) nin 2 . 



(15) 



In what follows we consider the case b > a. The Eq. 
([T5jl implies that rT\ — rvi = N is the integral of mo- 
tion for Eq. (|13p . Further we will consider N as a fixed 
integer(positive,zero or negative). Let us prove now that 



the dynamics of Eq. (|13|) can be considered as super- 
position of independent dynamics for each iV respec- 
tively . Indeed, let N is some positive integer and take as 
generating function the normalized function of the form: 
Gjv (u, v, t) = u N f (uv, t) (where / (1, t) = 1). It is easy 
to see that for such GF nT — rtii = TV for arbitrary 
f (x,t) = f (uv,t). Substituting this expression for GF 
in Eq. (|13p after simple algebra we obtain the closed 
equation of evolution for the function / (x, t): 



J 



df (x,t) 

at 



(i 



(N + l) — (b- ax) f + x— (b - ax) f 



(16) 



r 



Let us prove now that / (x, t) tends to b ^ ax when t tends hence satisfies to the equation: 
to infinity.For this purpose we introduce auxiliary func- 
tion h(x,t) = (b — ax) / (x,t). Then Eq. (fTB"]) implies 
the next equation for function h {x, t): 



dG dG 
u— v— = NG. 

dt dt 



(18) 



dh 
~dt 



(1 — x) (b — ax) 
(1-x) (b 



-JV 



dh 
dx 

ax) SF [h (x)} 



Sh 



d 2 h 
dx 2 



(17) 



The Eq. (|18|) implies relations which are similar to the 
relations (|11[) for foregoing model, namely: 



plies that 



<1F 
dt 



x)(b-ax) (^f) 2 dx < and 



where F [h (x)] is the functional of h (x) which has the 
i 9 

form: F[h(x)} = \ J x N+1 (||) dx; The Eq. (JT7J) im- 

o 

-/(I 



hence when t tends to infinity functional F [h (x, t)] tends 
to zero and therefore h (x, t) tends to certain constant 
value (remind that we assume b > a).QED.It is clear 
that all above arguments remain valid also in the case 
when integer N < 0, if one take as generating function 
G N = v\ N \g(uv,t) .Thus we reach the required conclu- 
siomthe general solution of Eq. p^|) admits the decom- 

N=oo 

position: G(u,v,t) — Kn Gn (u,v,t) . It means 

JV=-oo 

that for the study of evolution of arbitrary initial GF 
Go = G (u,v,t — 0) enough to consider separately dy- 
namics of each it component. Also we have proved that 
any component of initial GF of the form u N fo (uv) (such 
GF we again connect with certain pure state) transforms 

to the distribution - .Similar result is valid for 

(o— auv) 

the initial GF component of the form v M go (uv) .Thus 
qualitative description of evolution for arbitrary initial 
GF Go is completed. In conclusion of this part we present 
two notable results relating to evolution of pure states 
in this model. Resultl. Let Go = u N f (uv) .We have just 
proved that at any time G (it, v, t) has the same form and 



m - n 2 



71x712 



n 2 - n\n<i. 



N, 



(19) 



from which we are able to obtain expressions for varia- 
tions of statistical distributions of n\ and n 2l namely: 



nin 2 and 



a i = nj - (rii) = &2 — n 2 - (n 2 ) = nin 2 
hence correlation coefficient k between two populations 
1 and 2 is : k 12 = Kl " 2 ~ wm = l.Thus we come to the 
conclusion that all pure states in this model are maxi- 
mally correlated. 

Result 2. Let us again assume that initial pure state 



has GF of the form: Go 



fo (uv). Then mean oc- 



cupation numbers of two populations are: n\ 
^U=iand^=^U =1 . 



N 



When t tends to infinity then G tends to G 



st _ {b-a)u 



b—auv 

and occupation numbers in stationary state are: rTf = 
N + 53- and = Let ^ < 1 then initial num- 

ber n2 is very small, but in situation when ratio x = % 
tends to one, the number of individuals in both popula- 
tions increases with time without any bound. However, 
the difference between occupation numbers TTj" and 
remains constant as before. 
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V. THE FORMATION OF COALITION 
BETWEEN TWO POPULATIONS OR SOCIAL 
GROUPS IN THE STRUGGLE AGAINST 
COMMON RIVAL 

In this part we consider more complex model of inter- 
action between three distinct populations which is, in a 
sence, the superposition of two foregoing models. 

We will focus on the situation when two feeble groups 
of individuals enter into an alliance to resist together the 
common rival. Thus in this model on the one hand there is 
a cooperation between two feeble groups but on the other 
hand there is mutual antagonism in relation to common 
enemy. Unfortunately, I find it difficult to bring vivid ex- 
amples of such coalition in biology (although they are 
apparently exist). On the other hand it is clearly that in 
social, political and military conflicts there are the great 



number of forming such coalitions. However, I include 
this example in the present paper mainly to demonstrate 
how one can construct more and more complex models 
based on simple ones. 

Let us now turn to the explicit mathematical formu- 
lation of this model .As before we are starting from the 
Lindblad equation ([2]) but this time we choose as the 
set of Ri the next two operators: R\ = ^J^o-i a> 2 a- 3 and 

R 2 = \J\ a i a t a t ■ We believe that operator R\ describes 
events when more stronger population with index 1 wins 
in the conflict and respectively the coalition formed from 
two populations with indices 2 and 3 fail, while the oper- 
ator R 2 takes into account the contrary events. The Lind- 
blad equation for the diagonal elements of density matrix 
p„ „ „ in this case takes the form: 



dt 



m • n 2 ■ n 3 p ni _ _ - m • n 2 ■ n 3 p v 



+ b 



ni-n 2 - n 3 p— „ 2 



m-n 2 - n 3 p ri 



(20) 



( we again use in Eq. (f20| the notation: rii = rii + 1 and 
rii = rii - 1). 

As before we convert the system of difference- 
differential equations Eq. (|20|) to equivalent differen- 
tial equation for generation function : G (u, v, w, t) = 

S Pn 1 n 2 n 3 univn2wn3 ■ Required equation for GF in 

m ,n2 ,™ 

this model reads as: 

dG d 3 

^ = {u - vw) dML [{au - hvw)G] - (21) 



The Eq. (|21j) implies the next equations of motion for 
mean occupation numbers nT ,n2 , n^: 



-gj- = an 2 n 3 {ni + 1) - bn\ (n 2 + 1) (n 3 + I) 
dn 2 dn 3 



-an 2 n 3 (m + 1) + 6ni (n 2 + 1) (n 3 + 1) (22) 
at at 



One can see that Eq. (|22|) imply the existence of two inte- 
grals of motion, that can be written as:l) nT+na = N +p 
and 2) TTj" + = N.(we will consider further only the 
case when N and p are positive integers). The presence 
of these integrals allows one to give complete qualitative 
description of evolution (|2"TT) for any initial generating 
function Go (it, v, w, t — 0) .To realize this intention we 
firstly find the stationary solutions of Eq. (f2~Tj). Equat- 



ing its r.h.s to zero we obtain required result: 

G°*(u,v,w,) = A{U,V) + B(U, ™ )+C{V,W) , (23) 

au — bvw 

(where functions A (u, v) ,B (u,w) ,C (v, w) must be de- 
termined in such a way to obtain for G st some polynomial 
expansion.lt is easy to see that after relevant choice we 
can rewrite Eq. (1231) in the required polynomial form: 



G st {u,v,w) = [K (v) + L(w)} 



(au) + {au) 



1 (bvw) + ... (bvw) 



N 



(24) 
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(where polynomials K (v) , L (w) are determined by ini- 
tial GF ) . Now we will seek a complete set of pure states 
with self-closed dynamics, on which arbitrary generating 
function can be expanded. To this end let us consider the 



generating function in the form: G p = v p f (u,vw,t) . 
Substituting this expression in Eq. (f2~T|) after simple al- 
gebra we obtain the closed evolution equation for the 
function / (it, x,t) = / (u, vw, t), namely: 



df , , 
dt K ' du 



X dx* [{aU 



bx) f] + (p- 



i) Tx - bx ^ f\ 



(25) 



It is easy to see that the general solution of Eq. (|25|) can 
be represented as some linear superposition of homoge- 
neous polynomials of N degree in variables x and it, that 
is / (u, x, t) = Pn {u, x, t) 

N 

Besides it is clear, if we take the generation function 
in the form G p = w p g (it, vw) we come to the similar 
result. Now we are able to claim that the general solution 
of Eq. (|2ip can be expended on these two classes of pure 
states. We can write this decisive for the further study of 
this model result as: 



G (it, v, w, t) 



[v p Pn (u 1 vw)+w p Rn (u,vw)] 



(in Eq. (|26[) summation is over all integers N and p, and 
Pjv (it, x) , Rn (it, x) are homogeneous polynomials of N 
degree in it and x = vw depending on p as well) . 

The result obtained allows one to restrict the study of 
the Eq. ([2~T]) to the case of pure states evolution only. So, 
let us consider the pure state which continually has the 
form: G p n (u, v, w, t) = v p Rn (it, vw, t) (where Rn is 
homogeneous polynomial of N degree ).We specify here 
only two notable properties of such states in this model. 

1) Cooperating groups being in pure states are fully 
correlated, while conflicting groups are fully anticorre- 
lated.To prove this result note that pure states in ques- 
tion satisfy to the equations: 



dG 



pN 



dG 



0u 



■ w- 



dw 



NG 



(27) 



and 



dG 



pN 



dv 



dGpN 
W ~dw~ = P pN ' 



(28) 



The Eq. (I27| implies following relations:a) n\ + 713 = 
7V,and b) TTi+rii (ni — 1) + rT\n^ = Nni from which one 
can find the variance of first group population, namely: 



a \ = n\ — (ni) = l\ni — (ni) — 71177.3 = n\ 77,3 — 77,1773 
.Selfsame result one can obtain for the variance o 3 and 
hence correlation coefficient fci3 = " 1 " 3 ~" 1 " 3 between 

-1. Similar algebra results in 



groups 1 and 3 is equal to 
that fci2 = — 1, and &23 = l.QED.Let us consider now 
another important question: how the properties of sta- 
tionary pure state G s p l N depend on coefficients a and 
b and also on values of integers N and p. We intend to 
analyze this problem in all details elsewhere and in this 
paper, as an illustration, consider only single special case 
when integer N tends to infinity while integer p remains 
finite. So, let stationary GF is: 



Gst 
pN — 



(au) N + (au) N 1 (bvw) 



(bvwf . (29) 



We want to understand how in this stationary state ra- 
tio of mean occupation numbers in populations 1 and 3 
depend on the coefficients a and 6With the help of GF 
from Eq 
in the form 



the required ratio = can be represented 



ni 
"3 



.ab 



N-l 



Na N + (JV -l)a N ~ 1 b+ 

Nb N + (N - 1) b N ~ 1 a + a N - 1 b 



Nx N + (N - 1) x N ~ x + ...> 
N+(N - l)x+ ....^ N - 1 

I 



*-h \- Ln f N (*)] 

N - [Lnf (x) 



(30) 



x = f and f N (x) = 
Let N tends to infin- 



(where we use the notation: 

1 + *+.... + * N = ^r 1 ). 

ity Using the Eq. 

to (2=1) JV when k > 1, and M u nu, =r- 

V x ) 5 713 (1 — >c)N 

when k < 1. Thus in the case when x > 1 we obtain: 



(|30|) it is easy to see that attends 
tends to 



77i — > ^V, 772 — > P and 773 — > while in the case k < 1 
we get :tt7 -> j^,n^ -> N + (p - and 773" -> N.In 

other words if N 3> 1 and x > 1 the first population be- 
cames dominating, contrary in the case x < 1 coalition 
dominates. On the other hand when N ~ p the behavior 
of the model becames considerably more complex and 
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requires a separate study. 

Let us briefly summing up the main results of present 
paper. Based on ideas of QTOS we propose the consis- 
tent approach for statistical description of open classical 
systems with integer variables. We proved that for broad 
class of open systems possesing specific restrictions on 
the form of their interaction with environment the LME 
actually reduced to the Pauli master equation for diag- 



onal elements of density matrix. This fact gives one the 
reason to use the LME for quantitative study of different 
problems relating both to statistical physics and to var- 
ious "soft' sciences such as ecology, sociology, economics 
and so on. 

I would like to acknowledge L.A. Pastur for useful dis- 
cussions of the results of this paper. 
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